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Abstract 

In decays of the Z°, the dominant mechanism for the direct production of char- 
monium states is the decay of the Z° into a charm quark or antiquark followed by its 
fragmentation into the charmonium state. We calculate the fragmentation functions 
describing the splitting of charm quarks into S-wave charmonium states to leading 
order in the QCD coupling constant. Leading logarithms of Mz/m c are summed up 
using Altarelli-Parisi evolution equations. Our analytic result agrees with the complete 
leading order calculation of the rate for Z° — ► tpcc. We also use our fragmentation func- 
tions to calculate the production rate of heavy quarkonium states in W^, top quark, 
and Higgs decays. 



Introduction 

Among the rare decay modes of the Z° predicted by the Standard Gauge Theory 
are ones whose final states include charmonium. Of particular importance are the 3 Si char- 
monium states J /if) and if)', since their decays into lepton pairs provide easily identifiable 
experimental signatures. The dominant production mechanism for if> and ip' is the decay of B 
hadrons; in fact, this serves as a signature for B hadron production in Z° decay. The direct 
production of if> and if)' is therefore important in Z° decays as a background to B physics. It 
is also of interest in its own right, since it involves both short distance and long distance as- 
pects of quantum chromodynamics (QCD). The production of a charm quark and antiquark 
with small relative momentum in Z° decay is a short distance process with a characteristic 
length scale that can range from 1/Mz to as large as l/m c . The subsequent formation of 
a bound state from the cc pair is a long distance process involving all the complications 
of nonperturbative QCD. The methods of perturbative QCD can be used to calculate the 
production rates provided that it is possible to systematically separate the short distance 
effects from the long distance effects. 

Most previous work on charmonium production in Z° decay JT|, 0, |3| has focused on 
short distance processes in which the cc pair that form the if) is produced with a transverse 
separation of order 1/Mz- Long distance effects involved in the formation of the bound 
state are factored into the nonrelativistic radial wavefunction at the origin R(0). The best 
example of a short distance process is Z° — > ipgg, which has a branching fraction of about 
10~ 7 . This small branching fraction can be partly attributed to a factor of |i?(0)| 2 /(m c Mf ), 
which represents the probability for a cc pair that is produced in a region of size l/(m c M§) 
to form a bound state. This probability factor suppresses the branching fractions for short 
distance processes by m^/M|, so that they can be neglected in the limit Mz/m c — » oo. 

As pointed out by Kiihn and Schneider ||, the direct production of charmonium in 
Z° decay will be dominated not by short distance processes but by fragmentation processes. 
The fragmentation mechanism is the decay of the Z° into a final state that includes a high 
energy quark or gluon, followed by the splitting of that parton into the charmonium state 
plus other partons. In the fragmentation mechanism, the c and c that form the charmonium 
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state are produced with a separation of order l/m c . The probability that they form a bound 
state is proportional to |i?(0)| 2 /m? c . The branching ratio for such a process is therefore not 
suppressed by the factor m 2 c /M z associated with short distance processes. The fragmenta- 
tion of a parton is described by a fragmentation function D(z, //), which gives the probability 
for a parton with invariant mass less than // to split into the charmonium state with lon- 
gitudinal momentum fraction z. It was recently shown that the fragmentation functions 
for the splitting of partons into heavy quarkonium states can be calculated using perturba- 
tive QCD [§]. The fragmentation functions D g ^(z,/J l ) and D g ^, Vc (z, fi) that describe the 
splitting of gluons into S-wave quarkonium states were calculate to leading order in a s at 
the scale \i = m c . They were evolved to larger scales /i by using Altarelli-Parisi evolution 
equations, which sum up leading logarithms of n/m c . The production of ip in Z° decay from 
the splitting of virtual gluons has been considered by Hagiwara, Martin, and Stirling ||, 
but they did not organize the calculation in terms of fragmentation functions and were thus 
unable to sum up leading logarithms of Mz/m c . 

The production rate of ip via the process Z Q — > ipcc has been calculated by Barger, 
Cheung, and Keung ]7| with a rather surprising result: it has a branching fraction of about 
10~ 5 . This is almost two orders of magnitude larger than Z° — > ipgg, in spite of the fact 
that both rates are the same order in a s . An explanation for the relatively large branching 
fraction of Z° — > ipcc was provided in Ref. ||, where it was pointed out that this process 
includes a fragmentation contribution that is not suppressed by a factor of m 2 c /M z . This 
contribution can be factored into the rate for the Z° to decay into a cc pair multiplied by 
the probability for the c or c to fragment into ip. 

In this paper, we calculate the fragmentation functions D c ^(z,fi) and D c ^, Vc (z, /i) 
for a charm quark to split into an S-wave charmonium state. The fragmentation functions 
at the scale [i = m c are calculated to leading order in a s (m c ). Altarelli-Parisi equations are 
used to evolve them up to the scale [i = Mz/2 appropriate for Z° decay. Our simple analytic 
result for Z° — > ipcc agrees with the complete leading order calculation of Barger, Cheung, 
and Keung. We also use our fragmentation functions to calculate the direct production rates 
for i\) in decays and for T in top quark and Higgs decays. 
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Z° Decay via Fragmentation 

The fragmentation contribution to the inclusive decay rate of the Z° into charmonium 
is the term that survives in the limit Mz/m c — ► oo. The general form of the fragmentation 
contribution to the differential decay rate for the production of a if> of energy E is 

dT(Z° -> ip(E) + X) = V / dz df{Z° -> + X, /i) A_^(z, /i) , (1) 

where the sum is over partons of type % and z is the longitudinal momentum fraction of the 
if> relative to the parton. The physical interpretation of (jl]) is that a if> of energy E can be 
produced by first producing a parton % of larger energy E/z which subsequently splits into 
a if) carrying a fraction z of the parton energy. The expression for the differential decay 
rate has a factored form: all the dependence on the energy E is in the parton subprocess 
decay rate dT, while all the dependence on the charm quark mass m c is in the fragmentation 
function -Dj_>^. To maintain this factored form in spite of the logarithms of Mz/m c that 
arise in perturbation theory, a factorization scale \x must be introduced. The dependence 
on the arbitrary scale fi cancels between the two factors. Large logarithms of E/fj, in the 

subprocess decay rate T can be avoided by choosing \x on the order of E. Large logarithms 
of \ijm c then necessarily appear in the fragmentation functions Z?j_»^(z, /i), but they can be 
summed up by solving the evolution equations || 

r 1 dy 

Li-g-Di^(z,n) = Y,J — p i^j( z /y,v) Dj^(y,/i) , (2) 

where Pi-,j (x, y) is the Altarelli-Parisi function for the splitting of the parton of type % into 
a parton of type j with longitudinal momentum fraction x. For example, the c — > c splitting 
function for a charm quark with energy much greater than its mass is the usual splitting 
function for quarks: 

P_( I , rt = ^^ I _ 5 - +2 *(l-x)j. (3) 

The boundary condition on the evolution equation (0) is the initial fragmentation function 
Di^(z,m c ) at the scale m c . As shown in Ref. |J, it can be calculated perturbatively as a 
series in aAm r ). 



We can easily count the order in a s for the fragmentation contributions to ip pro- 
duction in Z° decay. The subprocess rate T for producing gluons is of order a s , while that 
for producing quarks is of order 1. The fragmentation function for a gluon to split into ip, 
which was calculated in Ref. [f|, is proportional to a^. A light quark can split into a ip 
only by radiating a gluon which splits into a ?/>, so its fragmentation function is of order . 
In contrast, the fragmentation function for a charm quark to split into a tp, which will be 
calculated explicitly below, is only of order a 2 s . Thus the fragmentation of charm quarks into 
ijj dominates by two powers of a s over the fragmentation of light quarks or gluons. 

Keeping only the charm quark and antiquark contributions to (p]), the energy distri- 
bution of the ip reduces at leading order in a s to 

dV 9 W 

— (Z° -> + X) = 2 T(Z° -> cc) (z, M z /2) ,z = — . (4) 

This fragmentation formula is of course applicable only for a ip of energy E that is a significant 
fraction z of the energy Mz/2 of the charm quark and much greater than the mass of the 
ip. In the factor of 2 accounts for the contribution from the fragmentation of the c. We 
have set the factorization scale /i to Mz/2 to avoid large logarithms from higher orders in 
perturbation theory. At leading order in a s , only the diagonal term in the evolution equation 
(H) survives: 

dy 
H^-D c ^(z,n) = / — P c ^ c (z/y,/i) D c ^(y,/i) . (5) 
Ofi Jz y 

Integrating (^) over the energy, the total rate for inclusive ^-production is 

T(Z° ^i) + X) = 2 F{Z° -> cc) ^ dz D c ^(z, m c ) . (6) 

We have set the fragmentation scale equal to m c by exploiting the fact that at leading order 
in a s the Altarelli-Parisi splitting function (§) satisfies Jq dxP c ^ c (x, /i) = 0. The evolution 
equation (^) then implies that the fragmentation probability Jq 1 dzD c ^(z, /i) does not evolve 
with the scale \i. 

Fragmentation function for c — > ip 
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We proceed to calculate the initial fragmentation function D c ^(z, m c ) for a charm 
quark to split into a ip to leading order in a s {m c ). Our strategy is to isolate the contribution 
Yi to the decay rate for Z° — > ipcc that arises from the fragmentation of the charm quark. 
We can then obtain the fragmentation probability Jq dzD(z) by dividing Ti by the rate T 
for Z° — > cc: 

To = / [«] (2tt)V(Z l^ol 2 , (7) 

where Z, q, and q are the 4-momenta of the Z°, c, and c, and [dq] = d 3 q/ (16n 3 qo) is 
the Lorentz- invariant phase space element. The square of the amplitude A Q for Z° — > cc, 
averaged over initial spins and summed over final spins and colors, is 



1 / Z a 7@\ ( 



m c ) , (8) 



where T a is the Z°cc vertex whose explicit form is not required. In the limit M z » m c , 
the factors of m c in the trace can be neglected. 

The rate for the decay Z° — > ipcc is 



Ti = 



2M, 



/ [dq] [dp] [dp'] (2vr) V(Z -q-p-p')\Y. \M* , (9) 



where q, p, and p' are the 4-momenta of the c, ip, and c. The four Feynman diagrams 
that contribute to the amplitude A\ at leading order in a s are shown in Figure 1. The 
contributions to the process Z° — > ipcc that correspond to the fragmentation of the charm 
quark come from the region of phase space in which the ip — c system has large momentum 
q = p+p' of order M z and small invariant mass s = q 2 of order m 2 c . To facilitate the extraction 
of the fragmentation probability, we write the 3-body phase space for the outgoing particles 
in an iterated form by introducing integrals over q and over s: 

[dq] [dp] [dp'] (2n) 4 5\Z -q-p-p) 

J ^f[dq][dq] {^) A 8\Z - q - q) J[dp][dp>] {2^5\q-p-p') . (10) 
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We also express the two-body phase space integral over p and p' in terms of the longitudinal 
momentum fraction z of the if). In a frame in which the virtual charm quark has the 4- 
momentum q = (q ,0,0,q 3 ), the longitudinal momentum fraction of the ip is z = {po + 
P3)/(Qo + 93) an d its transverse momentum is p± = (px,P2)- Expressed in terms of these 
variables, the Lorentz invariant phase space element is [dp] = dzd 2 p±/(167i 3 z). Integrating 
over the 4- momentum p' and over p±, the 2-body phase space integral reduces to 

/ [dp][dp'\ {2-nf5\q -p-p') = 1- fdzQL-^- -™L) . (11) 
J on Jo \ z I — z J 

We have set = 2m c , which is accurate up to relativistic corrections. If s = q 2 is of order 
m 2 , the delta function 5 4 (q — p — p') constrains p± to be of order m c . From the mass-shell 
condition, the component po — ^3 = (pi + 4m^)/(p +p 3 ) is of order m 2 /M z . Thus, to leading 
order in m c /Mz, we can set p = zq. 

We proceed to isolate the contribution to the amplitude A\ from the fragmentation 
of the charm quark. In covariant gauges, this contribution comes from both of the diagrams 
in Figure la and lb, while the diagrams in Figures lc and Id contain contributions from 
c fragmentation. In the axial gauge associated with the 4-vector q, the contribution from 
fragmentation of the charm quark comes only from the diagram shown in Figure la. The 
amplitude for Figure la in this gauge can be reduced to 

q-(2q-p) J 

We have used standard covariant Feynman rules JI[ for projecting the amplitude for produc- 
tion of a cc pair with equal 4-momenta p/2 onto the amplitude for production of a ip with 
4-momentum p. The parameter -R(O) is the value of the nonrelativistic radial wavefunction 
at the origin. Averaging over inital spins and summing over final spins and colors, the square 
of the amplitude reduces to 

U lll9 1287raMo)| 2 1 / - Z a Z f3 \ / ._, . \ 
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where D is a Dirac matrix that depends on q, q, and p. We need only keep the terms in D 
for which the Dirac trace in ( |i~3l) is of order m\M\. While q, q, and p all have components of 
order Mz, s = q 2 is of order m\ in the fragmentation region. Simplifying the Dirac matrix 
by dropping terms which are suppressed by powers of m c /M z , it reduces to 

D = (s 2 -2m 2 c s -47m 4 c ) 4 - (s -m 2 c )(s - 9m 2 c ) i> 

2 

+ 4 ^-y ((s + 7m 2 ) q - q p 1 — is — 5m 2 ) q ■ p p 1 — 8m 2 c q ■ q 



2 \ 2 
s — mz ■ 



+ 12 1^W^)J < 14 > 

We have exploited the fact that while p 1 and are both of order Mz, their product is 
only of order m c Mz- The coefficients of p 1 and ^ in (|TJ|) are all manifestly of order m 4 , so 
we can substitute p = zq for all the remaining factors of p. The Dirac trace in flI2|) is then 
proportional to tr(r a ^ Tptf)- It is now easy to divide I\ by the decay rate T given in (|7|) 
to obtain the fragmentation probability: 

8o2|i*(0)| 2 f» 1 / 4m 2 m 2 c \ 



z 1 — 



/ z(l-z) 

I (s 2 - 2m 2 s - 47m*) - z(s - m 2 c )(s - 9m 2 c ) + 4^ — -^s(s - m 2 c ) 

- 4 8 ~f_~/" 2 m^( g - m 2 c ) + 12^^(, - m c 2 ) 2 ) - (15) 

Note that the upper limit on the integral over s has been increased to oo. Since the integrand 
behaves like 1/s 2 at large s, this only changes the integral by an amount of order m 2 /Mf, 
which we have been systematically neglecting. Evaluating the integral over s in (|i"5[) , we 
obtain our final expression for the initial fragmentation function: 

64 2 \R(0)\ 2 z(l-z) 2 (16-32z + 72z 2 -32z* + 5z 4 ) 
D c ^(z,3m c ) = — a s (3m c ) — — . (16) 

We have set the scale \x in the fragmentation function and in the running coupling constant 
to fi = 3m c , which is the minimum value of the invariant mass ^/s of the fragmenting charm 



quark. We have also set 2m c — > My, in the denominator, which is accurate up to relativistic 
corrections. Integrating over z, we obtain the total fragmentation probability: 



Jl dz D c ^(z,3m c ) = ^^(3m c r l -^f- (H^-571og2) . (17) 



Fragmentation function for c — > rj c 

The fragmentation function for a charm quark to split into the 1 So state of charmo- 
nium f] c can be calculated in the same way as for tp. The starting point is the expression (|9]) 
for the decay rate for Z° — > r/ c cc, except that the amplitude A\ in (|T2[) must be replaced by 

A ' = « ^ «( z ) ^ '■( ! "*(73W S(^ '' ) ( W + 4mc)75W + mc, 

c — 777 ^ \ 

+ q-(2q-p) 175 ^ + 2rn c )jn(q) . (18) 
The square of the amplitude has the form (|T3|), except that the Dirac matrix D reduces to 
D = (s + 3ml)(s — hm 2 c ) 4 — (s — m 2 c )(s — 9m c ) p' 



2 

+ 4—— ^-r (( s - m c) Q ■ Q -(s-3m 2 c )q-p 

q-{2q-p) V 

— 777,^ \ ^ 

——-^-\q.pq.(q-p)^. (19) 

Following the same path as in the ip calculation, we find that the initial fragmentation 
function for r) c is 

Integrating over z, the fragmentation probability is 

Jl dz D c ^ c (z,3m c ) = ^^m c f^f - 371og2) . (21) 
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Decay of Z° into C harmonium 

From (P), the branching ratio for the decay of the Z° into ip relative to the decay into 



r Z ^ Vcc , ,o \R(0 2 , , 

= 0.1870 a s (3m c 2 1 22 

-> cc) v Ml K ' 

The value of the parameter -R(O) can be determined from the ip electronic width to be 
\R(0) | 2 = (0.82 GeV) 3 . Taking a s (3m c ) = 0.23, we find that the branching ratio (H) is 
1.8 x 10 -4 . The simple result ( p2"D agrees with the complete leading order calculation of 
Z° — > ipcc in Ref. after taking into account the differences in the values of R(0), a s , 
and the charm quark mass. A slightly larger value for the wavefunction at the origin was 
used in Ref. 0: |,R(0)| 2 = (0.92 GeV) 3 . It was also assumed implicitly in Ref. that 
Z° — > ipcc is a short distance process, so the running coupling constant was taken to be 
a s (Mz) ~ 0.15. As we have shown, the dominant contribution comes from a fragmentation 
process, and the appropriate scale of the coupling constant is definitely on the order of m c . 
Finally, instead of M 3 in the denominator of (0), the authors of Ref. [0] used (2m c ) 3 with 
m c = 1.35 GeV. The difference between and 2m c is a relativistic correction, which we 
have consistently ignored in this analysis. Corrections to the fragmentation approximation 
are of order (2M^/ Mz) 2 or about 0.4%, which is much smaller than the size of relativistic 
corrections and higher order perturbative corrections. 

The rate for production of r) c by fragmentation differs by less than 3% from that for 
ip. From (|2"ID , we obtain 

This agrees with the calculation of Ref. after taking into account the differences in the 
values of a s , R(0), and m c and an apparent algebraic error of a factor of 3. 

The energy distribution of the ^'s produced by the fragmentation of charm quarks 
in Z° decay is given in ([|). It is proportional to the fragmentation function evaluated at 
the scale Mz/2. The initial fragmentation function ([16]) at the scale 3m c is shown as a 
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solid line in Figure 1. It must be evolved up to the scale Mz/2 using the Altarelli-Parisi 
equation ([|) in order to sum up the leading logarithms of Mz / m c from higher order radiative 
corrections. The result is shown as the dotted line in Figure 2. The evolution softens 
the energy distribution, shifting the peak in the fragmentation function from z = 0.75 to 
z = 0.68. The energy distribution shown in Figure 2 should be accurate provided that 
the energy E of the ip is large compared to its mass, or equivalently z » 0.07. The 
fragmentation function for r] c production is also shown in Figure 2. It has a slightly softer 
distribution, but its behavior is otherwise similar to that for the ij). 

The expression ( p2"D also applies with minor modifications to the corresponding branch- 
ing ratio for T production: 

where -R(O) is the radial wavefunction at the origin for the T, which is determined from its 
electronic decay rate to be |i?(0)| 2 = (1.72 GeV) 3 . Taking a s (3m(,) = 0.17, we find that 
the branching ratio Q23| ) is 3.3 x 10 -5 . The fragmentation approximation for T production 
in Z° decay is not as accurate as it is for ip production. Corrections are on the order of 
{2M r /M z ) 2 , which is about 4%. 

Decay of W ± into ip 

About 1/3 of the decays of the W + will proceed through the channel W + — > cs. The 
mass of the W is sufficiently large that the dominant production mechanism for charmonium 
will be W + — > cs, followed by the fragmentation of the charm quark into charmonium. 
Fragmentation of the strange antiquark into ip is suppressed by a factor of a 2 s . The branching 
ratio for decay into ip relative to decay into cs is therefore smaller than (|22|) by a factor of 2: 

r (^ + - ^) - o 0935 a (3m ? W (25) 

Numerically this branching ratio is 9.2 x 10 -5 . Our analytic calculation of the fragmentation 
contribution is consistent with the full leading order calculation of Ref . . 

Decay of Top Quark into T 
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The top quark will probably decay almost exclusively into W + b. If the top quark 
is heavy enough, the dominant production mechanism for bottomonium in top quark decay 
will be t — > W + b, followed by the fragmentation of the b quark into bottomonium. The 
branching fraction for the direct decay into the 3 Si state T is one half of (f 



which has the numerical value 1.6 x 10~ 5 . The complete leading order calculation of the 
rate for t — > W + bT gives a branching fraction of 4 x 10~ 7 for a top quark with a mass of 
100 GeV J7|. The fragmentation formula (^) does not apply to such a small value of the 
top quark mass, since the maximum momentum of the T — b system is only 13 GeV, too 
small for the decay rate to be dominated by fragmentation. The simple result fl26|) is a good 
approximation if the mass of the top quark is closer to 150 GeV. 

Decay of Higgs into T 

If the Higgs mass is below the threshold for decay into W pairs, than its dominant 
decay mode will be H — > bb. The dominant production method for bottomonium in Higgs 
decay will be H — > bb, followed by the fragmentation of the b quark or antiquark into 
bottomonium. The branching fraction for the direct decay into the T is twice (^), because 
both the b and b can fragment into T: 

F(H -> Tbb) . x9 IMO)! 2 , , 

J___) = „, 870a , ( 3^mi, (27) 

This branching ratio is 3.3 x 10 -5 , which is probably too small for this decay mode to be 
useful cLS cL si gnal for an intermediate mass Higgs boson. 

Conclusions 

We have shown in this paper that the dominant mechanism for the direct production 
of charmonium in Z° decay is fragmentation, the production of a high energy charm quark or 
antiquark followed by its splitting into the charmonium state. Most previous calculations of 
charmonium production have considered only short-distance production mechanisms which 
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are suppressed by a factor of m^/M|. We calculated the fragmentation functions D(z, /z) for 
charm quarks or antiquarks to split into S-wave charmonium states to leading order in a s . 
The fragmentation functions satisfy Altarelli-Parisi evolution equations which can be used to 
sum up large logarithms of Mz/m c . These fragmentation functions are universal, applying 
to the production of heavy quarkonium in any high energy process that can produce heavy 
quarks with energy large compared to their mass. We applied them to the production of 
charmonium and bottomonium in decays of the Z°, W^, top quark, and Higgs boson. 

A complete calculation of the rate for ip production in Z° decay must include the 
production of the P-wave charmonium states Xcj, followed by their radiative decays into ip. In 
calculating the production of P-wave charmonium states, there are two distinct contributions 
that must be included at leading order in a s . The P-wave state can arise either from the 
production of a collinear cc pair in a color-singlet P-wave state, or from the production of 
a collinear cc pair in a color-octet S-wave state ||. The calculation of the fragmentation 
functions for the splitting of charm quarks into P-wave charmonium states will be presented 
elsewhere fpLOjl . 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant DE-FG02-91-ER40684. 
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Figure Captions 

1. The four Feynman diagrams for Z Q — > ipcc at leading order in a s . 

2. The fragmentation functions D c ^(z, ji) and D c ^ Vc (z, /i) as a function of z for fj, = 3m c 
(solid lines) and /i = Mz/2 (dotted lines). 



13 



